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Abstract 

We discuss recent results obtained for the Hamiltonian Mean Field model. The 
model describes a system of N fully-coupled particles in one dimension and shows 
a second-order phase transition from a clustered phase to a homogeneous one when 
the energy is increased. Strong chaos is found in correspondence to the critical point 
on top of a weak chaotic regime which characterizes the motion at low energies. For 
a small region around the critical point, we find anomalous (enhanced) diffusion and 
Levy walks in a transient temporal regime before the system relaxes to equilibrium. 
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1 Introduction 

In this contribution we present a novel study of chaos and anomalous diffusion 
in the dynamics of the single-particle motion in the Hamiltonian Mean Field 
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(HMF) model. This model describes N classical particles moving on the unit 
circle and interacting through an infinite range potential. It has been recently 
studied both analitycally and numerically [1-4]. At low energy the particles 
are clustered and, as the energy is increased, the system undergoes a second- 
order phase transition, reaching a high temperature "gaseous" phase where the 
particles are homogeneously distributed on the circle. The motion is chaotic 
and the Lyapunov exponent reaches a maximum at the critical point. 

Particle motion is superdiffusive for a range of energies close to the critical 
point, but only in a transient quasi-stationary regime preceeding equilibration. 
Particles can evaporate from the cluster showing Levy walks. Close to the crit- 
ical point, if the continuum (A^ — * oo) limit is performed before the t — > oo 
limit, canonical equilibrium is never reached and quasi-stationary superdif- 
fusive states live forever [1-4]. Though anomalous diffusion and Levy walks 
have been widely studied in the recent literature in models with a few degrees 
of freedom [5,6], there are very few studies for many-degrees-of- freedom sys- 
tems [7,8]. The HMF model has the advantage of having an exact solution 
in the canonical ensemble, and therefore microscopic dynamics can be put in 
connection to macroscopic features, thus offering a very interesting perspec- 
tive [3,4]. 

In the following sections we remind the reader the details of the model and 
then we discuss superdiffusion in connection to microscopic chaotic dynamics. 



2 The model 

HMF describes a system of classical particles moving on a circle, charac- 
terized by the angles 9i and the conjugate momenta pi. Each particle interact 
with all the others according to the following Hamiltonian: 

H{9,p)^K + V, (1) 

where 

i — 1 i,j — 1 

are the kinetic and potential energy. If we define a spin vector associated to 
each particle rrij = [cos{9i), sin{9i)] and a total magnetization M = -i: Y^iLi 
the Hamiltonian then describes A^ classical fully coupled spins, similarly to 
the XY model. With this interpretation, the system is a ferromagnet at low 
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Fig. 1. Total magnetization M as a function of the energy per particle U = E/N . 
Numerical simulations are compared with the canonical prediction (full curve). The 
insets show a schematic pictorial view of the condensed and the homogeneous phase. 
See text. 

energy and shows a second-order phase transition at the critical en- 
ergy Uc = Ec/N = 0.75 (corresponding to the critical temperature Tc = 0.5 
in the canonical ensemble) [1-3]). 

We show in fig.l the total magnetization vs. energy U = E/N. Numerical 
simulations (dots) are compared with the canonical ensemble result (full line) . 
Pictorial views of the low energy ferromagnetic (clustered) phase and high 
energy paramagnetic (gaseous) phases are shown in the insets. The equations 
of motion for the N particles are 

Oi=Pi, Pi = -Msin{9i - (f>) , i = l,...,N , (3) 

where (M, 0) are respectively the modulus and the phase of the total mag- 
netization vector M. These equations are formally equivalent to those of a 
perturbed pendulum. To study relaxation to canonical equilibrium, we solve 
these equations on the computer using fourth order symplectic algorithms (the 
details can be found in Refs [2,3]). We start the system in a given initial state 
and we compute 6i,Pi at each time step, and from them the total magneti- 
zation M and temperature T (through the relation T = 2 < K > /N). We 
consider systems with an increasing size N and different energies U = E/N. 
The single-particle motion is given by Eq.( 3), if we consider the mean-field 
pair (M , 0) as an external time dependent forcing. In the next section we focus 
our attention on the dynamics of the single-particle motion and we show that 
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Fig. 2. Single-particle phase space for three different energy regimes: U = 0.1 (< Uc), 
U = 0.69 (close to Uc) and [/ = !(> Uc)- See text for more details. 

it is chaotic and superdiffusive. 



3 Chaotic motion and superdiffusion 

In this section we discuss the numerical results. In fig. 2 we show single-particle 
phase spaces for three different energies: below (2a), close to (2b) and above 
(2c) the critical energy Uc = 0.75. While for U = 0.1 the particle never escapes 
the cluster and moves on the chaotic web near the center of the nonlinear 
resonance, close to the critical point (at Uc — 0.69) it escapes frequently from 
the cluster and is then trapped again, exhibiting the well known phenomenon 
of separatrix crossing. Above the critical energy one has mainly a ballistic 
regime, separatrix crossings are rarely observed and particles interact very 
weakly. 

We show in fig. 3 the largest Lyapunov exponent A^p for the single-particle 
motion as a function of the total energy per particle. The motion is chaotic 
and has a peak at the critical point, but Xgp is much smaller than the largest 
Lyapunov exponent calculated for the full dynamics in Refs. [2,3]. The reason 
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Fig. 3. (a) We show the largest Lyapunov exponent for the single-particle as a 
function of U. In the inset we show the behavior at small energy in a log-log plot, 
see text, (b) We report numerical simulations (open circles) that show that A^p 
scales like A^~^/^ (straight line). 

is that in the calculation of the Lyapunov of the full system one gets the addi- 
tional contribution to the instability of the off-diagonal terms of the Jacobian. 
We show in the inset a log-log plot with a fit at energies less than 0.3, which is 
a power law, i.e. A ~ U^^^; again quite far from the law found for the full 
dynamics. A very weak N dependence is present is this region. In the high en- 
ergy region A^p vanishes with A^, with the "universal" law A^p ~ A^~^/^, which 
is found in many different contexts [9]. This result is displayed in fig. 3(b) for 
U = 2. 

Diffusion and transport of a particle in a medium or in a fluid flow are char- 
acterized by the mean square displacement cr^(t) in the long-time limit. In 
general, one has 

a\t) ~ r (4) 



with a = 1 for normal diffusion. All the processes with a ^ 1 are termed 
anomalous diffusion, namely subdiffusion for < a < 1 and super diffusion for 
1 < a < 2. 

In order to study anomalous diffusion in HMF we follow the dynamics of each 
particle starting the system in a "water bag" , i.e. a far-off-equilibrium initial 
condition obtained by putting all the particles at 6*^ = and giving them a 
uniform distribution of momenta with a finite width centered around zero. We 
compute the variance of the one-particle angle 6 according to the expression 

a\t) =< {e- <d>f> , (5) 
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Fig. 4. (a) Levy walks for different test particles at U = 0.69. (b) Variance of the 
single-particle displacement as a function of time for four different energies. 

where < . > indicates the average over the N particles, and we fit the value 
of the exponent a in Eq. (4). 



In fig. 4(a) we show the motion of four particles in the same simulation as 
a function of time for = 500 and U = 0.69. The figure shows that the 
erratic motion is not brownian and Levy walks are present. Each particle 
experiences in time localized motion and free walks. In fig.4(b) we show the 
mean square displacement as a function of time. For energies close to the 
critical one the exponent a differs from 1, as it should be for normal diffusion, 
and takes values close to 1.4. One has therefore super diffusion, which however 
turns again into normal after a crossover time. This behavior does not seem to 
depend strongly on energy, although the super diffusion exponent a is slightly 
energy dependent. In Ref. [4] we have shown that this crossover time coincides 
with the equilibration one and that by using the model of Ref. [6] it is possible 
to interpret this dynamics in terms of Levy walks. 



4 Conclusion 



We have discussed chaotic and superdiffusive motion in a N-body Hamiltonian 
system which exhibits a second-order phase transition. The single-particle mo- 
tion is shown to be chaotic for energies smaller than the critical one. Superdif- 
fusion is present in a transient out-of-equilibrium regime for energies close 
to the critical one, where the system is strongly chaotic and quasi- stationary 
states exist. The crossover time from enhanced diffusion to normal coincides 
with the equilibration time which diverges with N. This is true for fully cou- 
pled systems, but should appear also if the interaction decays with the dis- 
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tance among the particles with a power law [10]. Our model lies in the class of 
self-consistent many-degrees-of-freedom dynamics, to which also belongs the 
single-wave model described in Ref. [11]. This is the closest contact of our 
results to the Levy- walk behavior found in fluid experiments [12]. 
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